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A nonextensive thermostatic approach to chaotic dynamical systems is developed by expressing 
generalized Tsallis distribution as escort distribution. We explicitly show the thermodynamic limit 
■ and also derive the Legendre Transform structure. As an application, bit variance is calculated for 

ergodic logistic map. Consistency of the formalism demands a relation between box size (e) and 
degree of nonextensivity, given as (1 — g) ~ — 1/ln e. This relation is numerically verified for the 
case of bit variance as well as using basic definition of Tsallis entropy. 



in 



o 



r-| , 05.45.+b, 05.20.-y, 05.70.Ce 

o 

I. INTRODUCTION 

tj^ ' Techniques borrowed from Gibbs-Boltzmann or extensive thermodynamics play important role in the characteri- 
zation of complex behaviour exhibited by dynamical systems. Analogous notions of entropy, temperature, pressure 
^ ■ and free energy can be applied to quantify the fractal or multifractal attractors of chaotic nonlinear mappings [||. 

In fact, thermodynamic theory of multifractals is applicable to any arbitrary probability distribution irrespective of 
' ' its mode of generation. A useful strategy here is the following: for the given probability distribution pi = exp(— 6^), 
where bi is bit number, form another set of normalized distributions Pi = Pi^ f^^Pi^, called escort distributions of 
order (3. For different /3 values, the new distributions are better able to scan different features of the original distribu- 
^ ' Hon Pi. Now in standard thermodynamics, maximizing the Shannon-Boltzmann entropy S — ~^^pi\ii pi under the 
constraint of given mean value of internal energy, one obtains the canonical distribution. Pi ~ exp(^' — /3Ei), where 
^ ^ = —In Z{l3) = —In iJ2i'^^Pi~(^^i)) is equivalent to Helmholtz free energy and Z{l3) is the partition function. 
^ In chaos theory, the escort distribution is also sought in a 'generalized' canonical form Pi — exp(\E' — /36i), where 
QQ ^I' ~ —In Z{f3) and Z{f3) = '^iPi^ = exp(— /36i). Thus in thermodynamic analogy, positive bit number bi is 
y—{ ' equivalent to energy Ei, (3 plays the role of inverse temperature and \1/ is a generalized free energy. Knowing the 
04 • partition function is the first step to the study of analogous thermodynamic quantities and relations in chaos theory. 
0\ ■ The analo gy goes deeper and interesting phenomena like phase transitions have been highlighted in dynamical systems 

On the other hand, it has been realized in recent years that the extensive formalism of thermodynamics fails to 
yield testable results (such as finite valued expressions for response functions of the various thermodynamic quantities) 
for systems with long-range interactions or which evolve with multifractal space-time constraints or have long term 
_ memory effects [Q. To deal with such cases, a nonextensive formalism of statistical thermodynamics was proposed by 
I . Tsallis This formalism has many successful applications to its credit by now, and the importance of nonextensive 
''O ' behaviour has been shown for diverse kinds of systems: Levy anomalous diffusion (6|j^, stellar polytropes 1^1, pure- 
C , electron plasma two dimensional turbulence [p|Jl0[|, solar neutrinos inverse brcmsstrahlung in plasma to 
O name a few [ p^ . Generally, this formalism is based on the proposal of a non- logarithmic entropy 

: s,{p) ^ ^-^f' ^ (1) 

^ Here, \bi] may just be considered as notation for the generalized bit number within the nonextensive framework. For 
q — > 1, \bi] bi — —In Pi and Eq. ([^) gives the Shannon entropy. Tsallis entropy is pseudo-additive (nonextensive) 
with respect to subsystems whose joint probabihties factorise: if pj^^^ — pjpV , then Sq{p^'^^^) = Sq{p^) -\- Sq{p^^) -\- 
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(1 — q)Sq{p^)Sq{p^^). Thus (1 — g) represents the deviation from extensivity. Maximising Tsalhs entropy under the 
internal energy constraint Uq — EiPi"^ yields the equilibrium Tsallis distribution 

^ {1 - (1 - g)/3gav(^-^) 

where f3 is the inverse temperature. For g — > 1, we get back the canonical distribution. 

It is believed that Tsallis formalism is a natural frame to study systems with fractal structures ||l5| . Indeed the form 
of Tsallis entropy was inspired from multifractal ideas |^ . Also its connection with a scale-invariant thermostatistics 
has been conjectured Simple relations exist between q and characteristic exponents of Levy flight distributions 
and anomalous diffusion IT^ . Recently, a precise relation was proposed between q and multifractal scaling properties 
of critical attractors of nonlinear one-dimensional maps JTst . For such cases, at the onset of chaos, the exponential 
sensitivity to initial conditions is replaced by power law sensitivity, and generalizaton of Kolmogorov- Sinai entropy 
has to be invoked on the lines of Tsallis formalism . 

Realizing the clear signatures of nonextensivity in chaotic dynamical systems, we feel a nonextensive thermodynamic 
approach needs to be developed for such systems. In this paper, we develop the thermostatistics of such a formalism. 
The plan of the paper is as follows: in section II, we present our formalism by first writing the partition function, 
then seeking its thermodynamic limit. By generalizing the Renyi dimensions, we establish the Legendre Transform 
structure. In section III, we present as an application, the calculation of bit variance within nonextensive framework. 
In section IV, we discuss the relation between box size (degree of coarse graining) and degree of nonextensivity (1-q). 
Section V is devoted to conclusions and outlook. 



II. GENERALIZED FORMALISM 



A. Basic ingredients 

Our approach is to work with escort distributions which are now sought not in canonical form, but as generalized 
Tsallis distribution. For this purpose, we consider a general, nonextensive statistics given by 

(3) 

For /3 — 1, this gives the relation between pi and generalized bit number [6,] in Tsallis statistics (see Eq. (|^)). Also for 
q — > 1, Pi = exp(— fei), which corresponds to extensive statistics. Now the escort distribution for the given distribution 
of Eq. (|) is 

p ^ _pI_ ^ {i-{i-q)mY'^'-'^ 
' n^p^' Eai-(i-9)/3[fo.]r/(i-^)' 

which is a generalized Tsallis distribution (see Eq. (2)), where now plays the role of energy Ei. The generalized 
partition function, Zq{f3) — J2iPi^ = ~ (1 ~ ■ Again for (3—1, escort distribution Pi is identical 

to original distribution pi. Note that in general, both the original distribution and the escort distribution depend on 
^ 

Now an important feature of Tsallis formalism is that it retains much of the formal structure of standard thermody- 
namics, including Legendre Transform structure Jl^ . This is facilitated by a generalization of the logarithm function 

as 

^^ = -\nqZq{(3)^-^^-l. (5) 

With the aid of such modified logarithm (or conversely, the exponential function), we develop the nonextensive 
thermostatistics of multifractals. Let us consider a certain nonlinear mapping over one-dimensional (for simplicity) 
phase space which has been partitioned into boxes of equal size e. Let pi represent visiting frequency of the map for 
a box i. Then we define local crowding indices as 

(1 - q)f3\n e 

which for g — > 1, can be written as ai (e) = IgiiL. From Eqs. (|) and (|), we may write 

[h]=a,ie,q)V, (7) 
where V — —In e, is the extensive volume parameter in thermodynamic analogy. 
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B. Thermodynamic limit 



To derive e — ^ or y ^ co limit of the generalized partition function, we replace the sum in Zq — J^iPi^ by an 
integral over a pO[|, 

^Qmax / \ 1/(1-9) 

Zg= da 7(a) 1 - (1 - q)f3aV , (8) 



where j{a)da are the number of boxes with in the interval [a, a + da]. Now in extensive case, 7(a) is assumed to 
scale as ^ = exp(/(Q!)F). In the nonextensive case, this can be generalized as 7(a) ~ (1 + {1 — q) f {a)Vy^^^^'^\ 

Thus we finally obtain 

Z,^ [ "^^^ da\ (l + il- q)fia)v] ( I - {1 - q)pav] 1 '\ (9) 



It is useful to remember that the intergrand above is in a sense, nonextensive generalization of the function exp((/(a) — 
f3a)V) in the extensive case. Now applying the saddle point method for the case of large V, we may write the above 
integral as 

Zg = ((l + il-q)fia)v)(l-il-q)f3av)\^'^ '\ (10) 



where a corresponds to the maximum value of the integrand for given P value. Then using Eq. (|^), we can write 

^ = {Pa- f{a)) - (1 - q)f3af{a)V. (11) 

Now in the limit V — > 00, the last term above diverges. To keep the generalized free energy density finite, we assume 
that 

(l-g)/3~l. (12) 



Then Eq. (|11| ) is given by 

lim^^{(3a-f{a))-af{a). (13) 

Note that the first term in parentheses is same as obtained for the extensive case. The second term arises due to the 
nonextensive nature of the present formalism. Also it may be noted that if we assume volume parameter V to be 
positive definite, then (1 — q)P > 0. 



C. Generalized Renyi Dimensions 

Next, we propose to generalize the usual Renyi dimensions [ pl| . This is required to discuss Legendre Transform 
structure of the present formalism. We define 

DM-I^It^I^JT.pA- (14) 



Again assuming the validity of Eq. (n2h for /3 7^ 0, in the e ^ limit, we get 



^.(/^) = Y^((E^^^')'"-i)- (15) 



For f3 — 0, the last term of the Eq. (|ll|) vanishes and thus, Eq. (|12| ) cannot be inferred. Then using Eq. (14) and 
V = —In e, we get Dq{0) = — lim^^o ^ , where N is the number of nonempty boxes. Also for /3 = 1, Eq. (|l|) gives 
Dq{l) = —(1 — q) J2iPd^Pi- To see the behaviour for large we write from Eqs. (||) and (0) 
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p/ = {l-(l-g)/3a,T/}i/(i~«). (16) 

Again in the limit of box size e — > 0, we assume the scaling form Eq. (|l^ ) and write 

p/ = {l-a,}i/(i-«). (17) 

For very large values of /3, the above expression is dominated by minimum value of a. Using this condition in Eq. 
(p^), we can write 

DM = ^«min, (18) 

which for (3 +oo yields Dq(oo) ~ amin- Similarly, for /3 — oo, one can show Dq{~oo) = amax- These results are 
thus identical to those obtained for the usual Renyi dimensions. 



D. Legendre Transform structure 

Now from Eq. ([Tsl), we write 



/3 

Using Eq. (|l^) with e = 0, and Eq. (|l9|), we get 

(/3 - l)i?,(/3) = /3a - (1 - a)/(a), (20) 

which we choose to write as 

T,(/3)=/3a-/,(a), (21) 
where we define Tg(/3) — {[3 — \)Dq{(i) and fq{a) = (1 — a)f[a). We find that 

%^ffl = ». (22) 

da op 

We note Tq{f3) being equivalent to the generalized free energy 'f, (that already obeys Legendre Transform structure), 
is a concave function of its argument. Then from Eq. (21), we infer that fq{a) is also a concave function, implying 
that the Legendre Transform structure is well defined. 



III. GENERALIZED BIT VARIANCE: AN EXAMPLE 



An alternative method to characterize multifractal statistics is using bit moments and bit cumulants of escort 
distributions p^ . Particularly, the second bit cumulant (r2) measures the variance of bit number. In the following, 
we derive expression for bit variance within the generalized statistical framework and apply it to ergodic logistic map. 
For nonextensive case, we write generalized second bit cumulant for escort distribution P as 

-(T.p^')'\T.^b^]p^'')y (23) 

The corresponding quantity for the original distribution (P — I) is given by 

^ ^=q[Y.ib.]V''''^ - (24) 



r,(.)(p) = -^ 
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Note that for g — > 1, we have T2{p) = Pi)'^ ~ iJ2iPi^'^ Pi)"^ ■ For the case of ergodic maps, bit variance density 

which is equivalent to heat capacity is given as C2 — ((In p)^) — (hi p}^ [ p3| , where p is the natural invariant density 
of the map. For the nonextensive case, we find 

C^^^-J^Mp^''-''^)-{p'?)- (25) 

In Fig. 1, we plot C2 and c'^^ against the nonlinearity parameter r for the logistic map, Xn+i — rxn{^ ~ Xn)- We 
find that for Q < q < 1, the bit variance (which represents size of fluctuations in bit number) c'^^ is larger than the 
quantity C2 of extensive case. We present an interpretation of this result in the next section. 

We also find quite an interesting feature i.e. the general trend of C2 vs. r curve, evaluated at smaller box size can 
be matched by C^^ evaluated at arbitrarily larger box size and appropriate value of 9 < 1. In Fig. 2, we show two 
cases to corroborate this point. Quite remarkably, the two quantities show matching for a wide range of the parameter 
r. To be precise, we note that decreasing the box size, lifts the heat capacity curve C2. Again, g < 1 also results in a 
general increase in heat capacity c'^^ (Fig. 1). Thus we can expect that heat capacity results with q = \ and smaller 
box size are reproduced with q < 1 and larger box size. This implies that a range of q and e values exist, which yield 
the same value for c'^^ . In the following, we make use of this possibility. 



IV. CONNECTION BETWEEN e AND (1 - Q) 

An important assumption in the above which makes the whole thermodynamic formalism consistent, is the relation 



between box size e and the degree of noncxtcnsivity (1 — q), as given by Eq. (12). Note that this relation is for an 
escort distribution of order /3. For the original distribution, (3=1 and we have (1 — g) ^ "1777- In the following, 
we numerically check the valdity of this scaling form. Again we consider the example of section III. For exactness, 
we now focus on a fixed value of nonlinearity parameter r = 3.81 in the chaotic region. The way we go about this 
exercise is as follows: for a given value of q and box size e, we calculate C2^'^\ Then we vary both q and e, so that 
we obtain same value of C2 within good approximation. Fig. 3 shows that a direct proportionality in fact exists 
between (1 — q) and — 1/ln e, for very small e. 

One can look for the origin of the proposed relation between (1 — q) and e, in the definition of Tsallis entropy 
itself. Thus consider a unit interval divided into W number of equal sized boxes. Assuming equiprobability, we have 

= ^ = e. Then Tsallis entropy for this ditribution is given by Sq = ^^^li ' ■ Again keeping Sq fixed, we vary q 
and e and as shown in Fig. 4, find that the proposed scaling form (1 — g) ^ 1/V is justified for large V . 

Finally, we make a remark on the results of Fig. 1. How do we interpret the higher values of (72^*^^ over C2 for g < 1 
in terms of effects of nonextensivity ? In the present context, we look at this issue from the viewpoint of inherent 
correlations due to nonextensivity. It is known that if the bit cumulants of subsystems, deviate from additivity, 
this indicates the presence of correlations. This happens naturally for the nonextensive case, where e.g. the first 
bit cumulant, F^*^ is identical to Tsallis entropy and hence is nonadditive with respect to statistically independent 
subsystems. Tsallis entropy is superextensive for < g < 1, i.e. the entropy of composite system is larger than the 
sum of entropies of individual subystems. This implies that negative correlations are set up for < g < 1. We propose 
that the overall increase in bit variance for the nonextensive case above, refiects these implicit correlations. In the 
thermodynamic limit of ^ — > 00, the proportionality between (1 — g) and \/V suggests that g ^ 1 — 0. Thus these 
correlations can be expected to be suppressed in this limit. Such correlations have also been studied in the context of 
classical ideal gas model within nonextensive framework [M . 



V. CONCLUSION AND OUTLOOK 



In the present paper, we presented a nonextensive generalization of thermostatistics of dynamical systems. We 
have seen that in the thermodynamic limit, generalized free energy density is finite if we assume direct proportionality 
between (1 — g) and 1/V. This relation also plays important role in exhibiting the Legendre Transform structure. 
We have numerically shown the validity of this relation for generalized bit variance and from the definition of Tsallis 
entropy itself. We know finite partitioning of the phase space is a fact of life in the study of chaotic systems. That 
this feature can be related to (1 — g), provides a new insight into the nature of nonextensivity. Further development 
of the formalism treating boxes with variable size, extension to dynamical apsects involving time evolution will be 
extremely welcome. It can be hoped that these possible developments will connect with the recent work ||l^, [fol on 
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low-dimensional dissipative systems and their power law sensitivity to initial conditions at bifurcation points or onset 
of chaos. 
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FIG. 1. Heat capacity equivalent to bit variance is plotted against control parameter r of logistic map. Solid curve represents 
C2 i.e. the extenive case, while dotted curve shows the behaviour of C^''^ . q value is set at 0.97. Higher value of bit variance 
in the latter case reflect the presence of nonextensive correlations. 
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FIG. 2. C2 (solid line) and (dashed line with symbols) can be matched for a wide range of r values: (a) C2 with 2^ 



boxes. Cj''' at q 



0.99 with 2^° boxes, (b) C2 with 2^^ boxes. Clf' at g = 0.98 with 2^" boxes. The plots show that decreasing 
q from unity in Cj'^ has the equivalent effect of decreasing the box size in C2. Also note that C2 increases in general, on 
decreasing the box size. 
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FIG. 3. Plot between (1 — q) and 1/V, keeping fixed C2 value originally calculated at r = 3.81, q = 0.9 and box size equal 
to 1/5000. The points can be fitted to straight line 1.1a; — 0.032. 
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FIG. 4. For equiprobablity, and pi = e, the plot between (1 — q) and 1/V, keeping fixed TsaUis entropy originally calculated 
at g = 0.9 and box size = 1/5000. Straight line fit is given by 2.26x — 0.17. 
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